Extreme wave events at the initial stage of the 
modulational instability in a forced and damped 
wave medium 

E V Belkin, A V Kirichok 1 , V M Kuklin 

Kharkov National University, Institute for High Technologies, 4 Svobody Sq., 
Kharkov 61077, Ukraine 

E-mail: ssLndyrcs@gmail.com 

Abstract. The modulational instability of waves in a medium under the action of an 
external monochromatic force and dissipation is considered. The model which describes 
the nonlinear stage of the modulation instability was constructed with using methods 
of so-called S'-theory. Following this theory, we take into account only interactions 
between modes with slowly varying relative phases. Within framework of this model, 
it is shown that the waves satisfying the criterion of abnormality appears as a rule 
at the initial stage of the instability. With further development of the process, the 
abnormal waves can be also observed but their amplitudes decrease gradually with 
decreasing of the average wave height. 
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1. Introduction 

In recent years, the attention of researchers was attracted to the phenomenon of large- 
amplitude short-lived waves observed in various nonlinear dispersive wave media. These 
waves was named as freaks or rogue waves. Of particular interest are the experimental 
observations of extreme ocean waves. A complete review on the various phenomena 
yielding to rogue waves in ocean can be found in the book pQ. Later, it has been 
found experimentally that freak waves can be generated in optical systems [2], [3] and 
in space plasma j3J. A possibility of generation of ion-acoustic freak waves in plasma 
with negative ions and generation of Alfven freak waves in astrophysical plasma was 
considered in [5]. 

Discussion about the possible mechanisms for the formation of freak waves has 
persisted for all last years. Now, the modulation instability is considered as main 
pretender on this role [6j. The phenomenon of modulational instability of finite 
amplitude wave was first investigated in [7], [S], [9]. In the conservative media, that 
is in the absence of external energy sources and mechanisms of energy absorption, 
the final stage in the development of modulational instability is the formation of a 
train of solitary waves or regimes with peaking (collapse) [10], [12J. That's why there 
are persistent attempts to describe the development of modulational instability using 
modified in various ways solutions relating to conservative systems (see, for example 
[13]). However, this approach may not necessarily be constructive, since it does not 
take into account the rather active and long transition process of the formation of 
stable wave structures. What is more, the stability of these structures can be achieved 
only in narrow ranges of parameters. 

Recently, Segur et al [H] and Wu et al [15] shown that dissipation stabilizes the 
modulational instability. In the wavenumber space, the region of instability shrinks 
as time increases. This means that any initially unstable mode of perturbation does 
not grow for ever. Damping can stop the growth of the side-bands before nonlinear 
interactions become important. Hence, when the perturbations are small initially, they 
cannot grow large enough for nonlinear resonant interaction between the carrier and the 
side-bands to become important. The amplitude of the side-bands can grow for a while 
and then oscillate in time. Kharif et al [16] considered the modulational instability of 
Stokes wavetrains suffering both effects of external forcing and dissipation. They found 
that the modulational instability depends on both frequency of the carrier wave and 
strength of the external force. 

This work shows that the presence of energy flow through such an open system 
results in complete or partial suppression or delay the process of forming a train of 
solitary waves and creates a new metastable state. This state, at least at the initial 
stage of the instability (i.e. during only a few times greater than the inverse linear 
increment of instability), characterized by the appearance of short-lived splashes of 
modulation of the fundamental wave with abnormally large amplitude. The average 
amplitude of the fundamental wave therewith remains large enough. In the case of a 
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large number of spatial modes in the excited spectrum, their amplitudes remain rather 
small in comparison with the amplitude of the fundamental wave and their interaction 
among themselves without the support of the fundamental wave can be neglected. At 
the initial stage of the instability, the width of the excited spectrum is rather wide 
and the process of interaction between side-band modes is accompanied in this case by 
so-called stimulated (i.e., governed by the fundamental wave) interference [17]. Further 
development of the instability at large times leads to a narrowing of the spectrum and the 
formation of stable soliton-like large-scale formations that resembles structures obtained 
in H3J. 

The method for describing the modulational instability in an open system was 
originally formulated in [18]. It was shown there that the pumping of spin waves by 
the homogeneous precession of magnetization (fc — > 0), that can be described by the 
Lighthill's equation, excites the symmetric spectrum satisfying the conditions of space- 
time synchronism of the form 2oj = oj(k) + u(—h) and 2k = = k — k. Growing 
side-band spectrum acts back on the pump that leads to "freezing" of its amplitude at 
the threshold level and restriction of the instability. 

Further improvement of the theory [TH] takes into account the weaker interaction 
between excited low-amplitude modes, and the main contribution to this interaction 
is provided by pairs of waves symmetric with respect to the pump u(k) + u(—h) = 
u(k')+u(—k'), which ensure that mentioned above conditions of space-time synchronizm 
are satisfied for all interacting modes. Later, Zaharov and co-authors (see detailed review 
[20] and book [IT]) have formulated the theory of the nonlinear stage of modulational 
instability based on this approach, which subsequently became known as S'-theory. 
Within framework of this theory they come to description in the language of correlation 
functions (AkA^,) = nxA(k — k') and {AkA^) = ax^(k + k'), taking into account 
the effect of total coupling (correlation) of phases <px for modes synchronously 

interacting with the uniform pump field and representing <jk = n.K-exp{— iipK}- 

Among important results of this extremely constructive S'-theory are the discovered 
dominant mechanism of instability saturation, which consist in back action of the excited 
spectrum on the pump wave under condition of low above-threshold level, and growth 
of influence of the phase mismatch on the saturation of the instability with increasing 
above-threshold level [21] . 

The model, considered below, is a continuation of the study [22], where we have 
analyzed the evolution of individual modes in the spectrum of plasma oscillations in the 
process of modulational instability. 

In present work we study more carefully not only the dynamics of individual 
modes but also the spatial-temporal evolution of wave packets (e.g. waves and their 
envelope). In addition, we analyze some specific futures of the instability, in particular, 
the symmetry breaking of the excited spectrum during the progress of the modulational 
instability in a medium with strong dispersion. This model, as shown below, allows to 
detect the formation of abnormally high peaks of field modulation and the appearance 
of extremely high waves. 
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2. Mathematical model 



Consider, as an example, the modulational instability of externally driven wave in a 
medium with sufficiently strong dispersion and weak dissipation, which can be observed 
in plasma wave-guides, as well as on the water surface and other physical situations. 
Weak dissipation, as noted in [Hj, provides necessary stabilizing effect on the modulation 
of fundamental waves, as an external source provides the permanent energy transfer to 
the wave motion. 

It is significant that here, such as in the above discussed cases, the amplitude of 
externally forcing wave retains large during the instability development and greatly 
exceeds the amplitudes of the excited spatial spectrum. Therefore, in conditions of 
finite level of dissipation the side-band modes excited by the modulational instability 
will interact with one other only if this interaction is supported by the fundamental 
wave . In other words, the dominant regime of interaction, at least in the early stage of 
the nonlinear regime of the instability is the interaction of modes symmetrically located 
in /c-space of the spectrum relative to the central mode of large-amplitude fundamental 
wave. 

We use below the following dispersion relation that characteristic, as an example, 
of gravity waves on deep water [23] : 

A 2 k 2 



gk<i + 



(i) 



where g is some dimensional coefficient (for the gravity waves on deep water it is the 
acceleration of gravity). 

The equations for the field amplitude which satisfies the dispersion equation Q 
can be represented as follows: 
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where A = 2ir/ko is the wave-length of the fundamental wave, 8 is the linear dissipation 
coefficient. 

Excluding the fundamental frequency Uq = y/gk : 
A K oc exp{— iu t + i(k + K)x + ifx} ~^A K oc exp{+i(k + K)x + ifx}, (3) 
we can rewrite Eq.Q: 
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Introducing real amplitudes and phases Ak = \uk\ exp(i<£>^) we have obtained the 
system of equations which describes the modulational instability in a medium with 
strong dispersion: 



du 



K 



dr 
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U- 



Kulsm^ K + u^ K 22 u K U- K sin($ K -$ K ) 



(4) 



where §k = 2<^o — — f-K is the total phase (or the phase of the instability channel). 
A distinction needs to be drawn between modes with wave numbers K and k and phases 
$x and $ K . 
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where the mode frequencies are u)(K) — uio = u)[ko(l + K)\ — co(k ) = gk (l + K) 
y/gk . We use the following notations &; £ = (, u t/2 = r/a, a = &q|A | 2 
tyWo{k \A \) 2 /2, K = (k- k )/k , u tK 2 /A = (r/2a)K 2 , A K /A (r = 0) = 
uk exp{^#}, and also 

A K = 2 (y(l + K) - 1] + y/{l-K) - l) /a, 

P K = 2(1 + K) 2b + 2(1 - K) 2 - 5 - 2. 

Note that radicands may be not expanded in further calculations. 

Equations for amplitude and phase of the fundamental wave take the form: 
du Q 



dr 

chpo 
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(7) 



Here G is an external source, supporting the monochromatic fundamental wave. The 
summing over K can be replaced by summing over K m = mAK/ko, where AK = 
2Kma,x/N and m = ±(1, 2, N) with K maiX /k = y/2k \A \ = V2a, N is the number of 
modes. Note that we consider direct external forcing rather than parametric as in [16]. 

Equations (|4])-([7|) form a closed system that should be solved numerically. 

The wave profile in £-space looks as follows: 
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Figure 1. Time evolution of wave height parameters: Ucp - the average height of all 
waves in the observation domain, Uswh - the average height of a highest third, Umax 
- the maximum wave height [25] . 



Here, the fundamental wave is at rest and modulation of moves. 

Since the dispersion is sufficiently strong in the vicinity of the fundamental wave 
frequency, one can readily see from above equations the symmetry breaking in excitation 
of the Stokes and anti-Stokes part of the spectrum, which can significantly change the 
character of the modulation. 

3. Numerical results 

In order to analyze the wave height distribution (e.g. the distribution of vertical 
distances between the wave crest and the deepest trough preceding or following the 
crest ), we take a third of highest waves. Then we find the average height of all 
waves Ucp, average height of a highest third Uswh and the maximum wave height 
Umax in consideration domain (( C I = 27i/(AK/k ) = ixN/K m = nN/\/2a, where 
AK = 2K max /N, ( = k x). Calculations were performed for 600 modes in the spectrum. 
The ratio of dissipation level S to the maximum growth rate was chosen as 0.1 (e.g. 
5 = 0.1). In order to provide the unit amplitude of the fundamental wave at the initial 
stage of the instability we have also chosen the value of the external drive force as 
G = 5 = 0.1 

Note that the criterion which defines the extremely high waves 

U AG > 2U.SWH (9) 
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(a) (b) 

Figure 2. Wave profile (see ([8])) near extremely high waves for various realization of 
the process: (a) t cx 10; (b) r oc 30 



or something like this, should be used with caution because this criterion as usual is 
applied to statistics obtained on sufficiently large observation periods, but the highest 
amplitudes are observed at the initial stage of the instability. 

Nevertheless, much smaller wave amplitudes or heights also satisfy this criterion 
in the regime of developed instability, since a decrease in both medium and large wave 
amplitudes is observed with time. Let give two examples of realizations where some 
waves satisfy the criterion ^ (see Figj2]). 

Analysis of observations and numerical simulation show [1] that extreme wave 
events often occur in wave groups having the form of soliton-like structures. 

With decreasing dissipation level S the processes of energy exchange between the 
side spectrum and the fundamental wave is amplified, which is clear seen in Figure [3j 

It can be seen that the maximum of instability growth rate shifts to lower values 
of K with decreasing amplitude of the fundamental wave at the initial stage of 
instability. Since the maximum growth rate corresponds to a value of the total phase 
<& K = 2ip — {p K — {p_ K equal to 7r/2, the majority of mode pairs with different values of K 
are synchronized with the phase of the fundamental wave [26J during instability progress. 
The fact that the total phase is focused near the ir/2 sets the stage for almost one- 
type interaction of different mode pairs with the fundamental wave. This collective 
interaction of the growing spectrum modes with the fundamental wave and explains the 
nature of intense energy exchange between the fundamental wave and the spectrum 
at the initial stage of instability. Later, the spread of total phases grows and the 
energy exchange between the fundamental wave and the spectrum is attenuated. Owing 
to the fact that the phases of individual modes tpx remain randomly distributed (in 
particular, no symmetry of the form ipx ^ <f-K, <-Pk ^ —<P-k), the instability spectrum 
synchronized with the fundamental wave forms a different interference structure for 
each realization. Nevertheless, the degree of coherence of excited modes achieves 
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Figure 3. The total intensity of side-band spectrum for different values of 5 (N = 100) 



the maximum at the initial stage of the instability. Incidentally, the phenomenon 
of synchronization of total phases in the spectrum of modulational instability in the 
presence of the symmetry <p K = ip-x, usually leads to regimes with peaking [27]. The 
absence of the phase symmetry in pairs of interacting waves could weaken the intensity 
of interference peaks and reduce the time of their existence. 

Considering the dynamics of instability spectrum, it is possible to detect the 
phenomenon of its noticeable shift in relation to the spectral domain of the linear 
instability. This shift can be explained by reduced amplitude of the fundamental wave. 
In addition, note that the amplitudes of individual modes of the spectrum remain much 
smaller than the amplitude of the fundamental wave. It should be noted also the 
asymmetry of the instability spectrum relative to the fundamental wave due to strong 
dispersion and large enough modulational instability growth rate for fundamental waves. 

Analyzing the instability spectra one can see that the modulation length increases 
almost 3.5 times with development of the instability. Evolution of the relative 
characteristic modulation length is shown in Figure [5] 

Dynamics of two-dimensional wave processes turns out to be similar. Thus, the 
number of waves on the modulation length at the initial stage of the instability is much 
less than at the later stages of the process. 

Analyzing the probability of extreme wave events for different realizations of the 
process, we have found that one such wave appears among the 10000 waves, which is 
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Figure 4. The spectrum of instability at r cx 10 (a), r oc 35 (b), t« 140 (c) 



qualitatively consistent with the known observations of ocean waves. 



4. Conclusion 



The considered model shows that the appearance of waves with abnormally large 
amplitudes is characteristic of the initial stage of the modulational instability. The 
average and maximum wave heights noticeably decrease with the development of the 
instability. However, even in the later case (r oc 30) abnormally high waves can be 
detected according to the criterion d9j), although their amplitude is already one and a 
half - two times less than in the most interesting case (r oc 10). 

There are few waves fits at the modulation length in the initial stage of the process 
and one of which may be abnormally large in some realizations (our attention to this 
fact was called by E.Pelinovsky, pQ). At the stage of developed instability, the number 
of waves for the modulation length increases to three or four times. The probability 
of abnormally high wave events in our numerical simulations is qualitatively consistent 
with the known observations of ocean waves. It should be noted that the total energy 
of wave motion is approximately conserved with reasonable accuracy in the domain 
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Figure 5. Evolution of the relative characteristic modulation length with development 
of the instability. 



2 000- 
1 500- 



1 000 ■ 




Figure 6. Distribution of wave amplitudes in different realizations at the moment of 
r oc 10. 

under consideration (both in time and space) and the condition that the amplitudes of 
individual side-band modes remain much smaller of the fundamental wave amplitude 
is satisfied during the entire simulation time. The latter allows us to assume such a 
description of the modulational instability sufficiently correct. 

The authors are grateful to Prof. E.A. Kuznetsov for his interest and helpful 
comments. 
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